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Abstract 



The Wigner phase-space distribution function provides the basis for Moyal's deformation quantization 
alternative to the more conventional Hilbert space and path integral quantizations. General features of 
time-independent Wigner functions are explored here, including the functional ("star") eigenvalue equa- 
tions they satisfy; their projective orthogonality spectral properties; their Darboux ( "supersymmetric" ) 
isospectral potential recursions; and their canonical transformations. These features are illustrated ex- 
plicitly through simple solvable potentials: the harmonic oscillator, the linear potential, the Poschl- Teller 
potential, and the Liouville potential. 



1 Introduction 

Wigner functions have been receiving increasing attention in quantum optics, dynamical systems, and 
the algebraic structures of M-theory They were invented by Wigner and Szilard [0], and serve as a 
phase-space distribution alternative to the density matrix, to whose matrix elements they are related by 
Fourier transformation. The diagonal, hence real, time-independent pure-state Wigner function f(x,p) 
corresponding to the eigenfunction tp of Hip = Eip, is 

f(x,p) = ^JdyiP*(x- \y) e-WfPix + \y). (1) 

These functions are not quite probability distribution functions, as they are not necessarily positive — 
illustrated below. However, upon integration over p or x, they yield bona-fide positive probability distri- 
butions, in x or p respectively. 

Wigner functions underlie Moyal's formulation of Quantum Mechanics || , through the unique [f|, || one- 
parameter (h) associative deformation of the Poisson-Bracket structure of classical mechanics. Expectation 
values can be computed on the basis of phase-space c-number functions: given an operator A(x, p), the 
corresponding phase-space function A(x,p) obtained by p i— > p, xhi yields that operator's expectation 
value through 

(A) = J dxdp f(x,p)A(x,p) , (2) 

assuming the usual normalization, J dxdpf(x,p) = 1, and further assuming Weyl ordering, as addressed 
by Moyal, who took matrix elements of all such operators: 

A(x, p) = f drdadxdp A(x,p) exp(ir(p — p) + ia(x — x)) . (3) 
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Wigner functions are c-numbers, but they compose with each other nonlocally. Properties of these 
compositions were explored in, e.g., || [jj], and were codified in an elegant system in ||: to parallel operator 
multiplication, the Wigner functions compose with each other through the associative star product 

+ = e ^(dxd P -dpdx) _ ^ 

Recalling the action of a translation operator exp(ad x ) h(x) = h(x + a), it is evident that the ^-product 
induces simple "Bopp" -shifts: 

f(x,p)*g(x,p) = f(x,p- y d x ) g(x,P+Y d x ) = f{x + y d p ,p - y d x ) g(x,p), (5) 

etc., where d and d here act on the arguments of / and g, respectively. This intricate convolution samples 
the Wigner function over the entire phase space, and thus provides an alternative to operator multiplication 
in Hilbert space. 

Antisymmetrizing and symmetrizing the star product, yields the Moyal (Sine) Bracket ||, 

{{f,9}} = ^ ' ^ 

and Baker's || Cosine Bracket 

{if, 9)) = g ' 

respectively. Note |8| that 

J dpdx f*g=f dpalx fg . (8) 
Further note the Wigner distribution has a *-factorizable integrand, 

f(x, -2p) = ^Jdy (fV) e WP ) * fa) e^) . (9) 



In general, systematic specification of time-dependent Wigner functions is predicated on the eigenvalue 
spectrum of the time-independent problem. For pure-state static distributions, Wigner and, more explicitly, 
Moyal, showed that 

{{H(x,p),f(x,p)}}=0, (10) 

i.e. H and / ^-commute. However, there is a more powerful functional equation, the "star-genvalue" 
equation, which holds for the time-independent pure state Wigner functions (Lemma 1), and amounts to 
a complete characterization of them (Lemma 2). 

We will explore the features of this *-genvalue equation, and illustrate its utility on a number of 
solvable potentials, including both the harmonic oscillator and the linear one. The ^-multiplications of 
Wigner functions will be seen to parallel Hilbert space operations in marked detail. The Poschl- Teller 
potential will reveal how the hierarchy of factorizable Hamiltonians familiar from supersymmetric quantum 
mechanics finds its full analog in *-space. We determine the Wigner function's transformation properties 
under (phase-space volume-preserving) canonical transformations, which we finally elaborate in the context 
of the Liouville potential. 



2 The *-genvalue Equation 



Lemma 1: Static, pure-state Wigner functions obey the -k-genvalue equation, 

H(x,p)* f(x,p) = E f(x,p) . 



(11) 
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Without essential loss of generality, consider H(x,p) = p 2 /2m + V{x). 

H(x,p)*f{x,p) = ^({p-i^ x ) 2 /2m + V(x)^ jdye-^ + ^ a ^r(x-\y) *P(x + ~y) (12) 
= ^J d y((P- 4 9x) 2 /^m + V(x + ^y)) e~ iyp ^*{x - ^y) ^(x + ^y) 

= h S dy e ~ WP ( (i ^ + 4 ^) 2 /2m + V{x + r(x - \y) Hx + \y) 

= ^J d y e- typ r(x - \y) E ip{x + ^y) = E f(x,p), 
since the action of the effective differential operators on ip* turns out to be null; and, likewise, 

f*H = ±-Jdy e-™ (-(d v ~\ d x ) 2 /2m + V(x - |y)) ^(x - \y) ^{x + ~y) = Ef(x,p). (13) 

Thus, both of the above relations ( |To|) and Lemma 1 obtain. ■ 

This time-independent equation was introduced in ref H, such that the expectation of the energy 
H(x,p) in a pure state time-independent Wigner function f(x,p) is given by 

J H(x,p)f(x,p) dxdp = E J f(x,p) dxdp. (14) 

On account of the integration property of the star product, (^), the left hand side of this amounts to 
Jdxdp H(x,p) * f(x,p). Implicitly, this equation could have been inferred from the Bloch equation of 
the temperature- and time-dependent Wigner function, in the early work of |J. *-genvalue equations are 
discussed in some depth in the second of refs M, and in [10]. 

By virtue of this equation, Fairlie also derived the general ^-orthogonality and spectral projection prop- 
erties of static Wigner functions [M]. His results were later formalized in the spectral theory of the second 
of refs H (e.g. eqn 4.4). Consider g corresponding to the (normalized) eigenfunction ip g corresponding to 
energy E g . By Lemma 1, and the associativity of the ^-product, 

f*H*g = E f f*g = E g f*g. (15) 

Then, if E g ^ Ef, this is only satisfied by 

/*5 = 0. (16) 
N.B. the integrated version is familiar from Wigner's paper, 

J dxdp f * g = J dxdp fg = 0, (17) 

and demonstrates that all overlapping Wigner functions cannot be everywhere positive. The unintegrated 
relation introduced by Fairlie appears local, but is, of course, highly nonlocal, by virtue of the convolving 
action of the ^-product. 

Precluding degeneracy, for f = g, 

f*H*f = E f f*f = H*f*f, (18) 
which leads, by virtue of associativity, to the normalization relation || , 

/*/«/. (19) 
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Both relations (|T^) and (|T^) can be checked directly: 

1T1 — * zfi * — 

f(x,p)*g(x,p) = f(x,p- — d x ) g(x,P + y d x ) (20) 

= WfI dydY e ~ t{y+Y)p r f {x ~\ y + \ Y) ^ f{x + \ y+ \ Y) ~\ Y ~ \ y) + \ Y ~ \ y) = 
(sf I d{Y + y) e ~ l{y+Y)P ^ x ~\( Y+ y ^ ^ x + \^ y+ y ») (jtj d( h{X 2 V) )r A {Y " y^A {Y " y)) ) ■ 

The second integral factor is or 1/h, depending on / ^ g or / = g, respectively, specifying the normal- 
ization / * / = f jh in (|l~9l). In conclusion, 

Corollary 1: fa*fb = \ $a,b fa ■ 

These spectral properties are summoned up by their own necessity; much of their meaning, nevertheless, 
resides in their margins: For non-normalizable wavefunctions, the above second integral factor may diverge, 
as illustrated below for the linear potential, but the orthogonality properties still hold. 

Thus, e.g., for an arbitrary function(al) F(z), 

F[f*]f = F(l/h)f, (21) 

and, for *-genfunctions of Lemma 1, 

F[H*]f = F(E)f. (22) 

Baker's converse construction extends to a full converse of Lemma 1, namely 

Lemma 2: Real solutions of H(x,p) * f(x,p) = E f(x,p) (= f(x,p) * H(x,p)) must be of the Wigner 
form, f = Jdy e~ typ ip*(x — ^y)ip(x + |y)/27T, such that Hip = Eip. 

As seen above, the pair of ^-eigenvalue equations dictate, for f(x,p) = fdy e~ wp f(x, y), 

I dy e ~ iyp {~L^ v 4 + v{x ± \ y) - E ) = °- (23) 

This constrains f(x,y) to consist of bilinears ip*(x — %y) ip(x + \y) of unnormalized eigenfunctions ip(x) 
corresponding to the same eigenvalue E in the Schrodinger equation with potential V. ■ 

These two Lemmata then amount to the statement that, for real functions f(x,p), the Wigner form is 
equivalent to compliance with the -k-genvalue equation (real and imaginary part). 



3 Example: The Simple Harmonic Oscillator 

The eigenvalue equation of Lemma 1 may be solved directly to produce the Wigner functions for specific 
potentials, without first solving the corresponding Schrodinger problem (as in, e.g. Following for 

the harmonic oscillator, H = (p 2 + x 2 )/2 (with h = 1, m = 1), the resulting equation is 

((* + \d P ? + (P~ \d x ) 2 - 2E^j f(x,p) = 0. (24) 

By virtue of its imaginary part, (xd p — pd x )f = 0, / is seen to depend on only one variable, z = 4i? = 
2(x 2 +p 2 ), so the equation reduces to a simple ODE, 

(^- z d 2 -d z -E^f(z) = 0. (25) 
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Moreover, setting f(z) = exp(—z/2)L(z), this yields 

zd 2 z + (1 - z)d z +E-~) Hz) = 0, (26) 

which is the equation satisfied by Laguerre polynomials, L n = e z d n (e~ z z n ), for n = E — 1/2 = 0, 1, 2, 
so that the un-normalized eigen-Wigner-functions are 

/„ = e~ 2H L n (4H); L = 1, L 1 = l-4H, L 2 = 16H 2 - 16H + 2, ... (27) 

Note the eigenfunctions are not positive definite, and are the only ones satisfying the boundary conditions, 
/(0) finite, and f(z) — > 0, as z — > oo. 

In fact, Dirac's Hamiltonian factorization method for algebraic solution carries through (cf. |J) intact 
in *-space. Indeed, 

ff = i(x-tp)*(x + ip) + -, (28) 



motivating definition of 

Thus, noting 
and also that, by above, 



1 1 

a=—^(x + ip), a) = -j={x - ip). (29) 

a* a) — a) *a = 1, (30) 

a*ib = ^(x + ip)*e- (x2+p2) =0, (31) 
v2 

provides a *-Fock vacuum, it is evident that associativity of the ^-product permits the entire ladder 
spectrum generation to go through as usual. The *-genstates of the Hamiltonian, s.t. H * / = / * H, are 
thus 

/„ oc (ah) n f (*a) n . (32) 

These states are real, like the Gaussian ground state, and are thus left-right symmetric *-genstates. They 
are also transparently ^-orthogonal for different eigenvalues; and they project to themselves, as they should, 
since the Gaussian ground state does, /o * /o ex /q. It will be seen below that even the generalization of 
this factorization method for isospectral potential pairs goes through without difficulty. 



4 Further Example: The Linear Potential 

For simplicity, take m = 1/2, fi = 1. Recall ( fi2[| ) that the problem readily reduces to a free particle: 
H(x,p) = p 2 + x i — ^ Hfree = -P is accomplished by canonically transforming through the generating 
function F(x, X) = -\ X 3 - xX. The energy eigenfunctions are Airy functions, 

1 r+oc 

ifj E ( x ) = — dX e iFi - x ' X) e iEX = Ai(x - E) . (33) 

2-7T J-oo 

The *-genvalue equation in this case is 

(V + \d p ) + {p- % -d x ) 2 - E^j f(x,p) = , (34) 

whose imaginary part, \\d p — pd x ^j f(x,p) = 0, gives f(x,p) = f(x + p 2 ) = f(H). The real part of the 
equation is then an ordinary second order equation, just as in the above harmonic oscillator case. Moreover, 
here the real part of the *-genvalue equation is essentially the same as the usual energy eigenvalue equation: 

(z-\# t -E)f{z) = Q, (35) 
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where z = x +p 2 . Hence, the Wigner function is again an Airy function, like the above wavefunctions, 
except that the argument has a different scale and shiftu. 

2 2 /3 2 2 / 3 1 r 

f(x, P ) = —Ai(2 2 l\z - E)) = —Ai(2 2 / 3 (x + p 2 - E)) = — j dy e^ E ~^ ~* ^ . (36) 

The Airy functions are not square integrable, so that the conventional normalization /*/ = f does not 
strictly apply. On the other hand, the energy eigenfunctions are non-degenerate, and the general Corollary 
1 projection relations f a * ft> oc S a ^f a still hold for the continuous spectrum: 

fEl*fE2 = fEl((x + ^d p )+p 2 )fE2({x- t -d P )+ P 2 ) 

_ . 1 J dy d y e iy(El-x-(p-Y/2) 2 -y 2 /12) e iY(E2~x-(p+y/2) 2 -Y 2 /12) 

= 1 Jd(y + Y) e «v+Y)(^-*-P 2 -(y+rr/u) Jd {y ~ Y) e^C*i-«0 
= J-6(E1-E2) fE^Eiix+p 2 ), (37) 
by virtue of the direct definition (B9). 



5 Darboux Construction of Wigner Function Recursions 

Analogous ladder operators for eigenstates corresponding to "essentially isospectral" pairs of partner po- 
tentials |l4[] (familiar from supersymmetric quantum mechanics) can also be defined mutatis-mutandis for 
Wigner functions and *-products. They faithfully parallel the differential equation structures. 

Consider a positive semidefinite Hamiltonian 

H = p 2 /2m + V(x). (38) 
This can be written as a ^-product of two operators, 

H = Q* -kQ = (-3= + iW(x)^j * (-J= - tWOr)) , (39) 

provided 

W 2 - -^=d x W = V{x) . (40) 
V2m 



This Riccati equation, familiar from ssQM, can be Darboux-transformed by changing variable for the 
"superpotential" W(x): 

W = - n *L*° , (41) 
which reduces the condition to the Schrodinger equation for zero eigenvalue: 

-lf-dli> o + V(x)i/j o = 0. (42) 
2m 

Also note Q * f = for the corresponding Wigner function. It is easy to generalize this by adding a 
constant to H to shift the ground state eigenvalue from zero. 



1 This case is similar to the Gaussian wave function, i.e. the harmonic oscillator ground state encountered above, whose 
Wigner function is also a Gaussian, but of different width. S Habib kindly informed us that this solution is also given in ref 
H], eqn (29). 
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By virtue of associativity, it is evident that the partner Hamiltonian 



H' = Q*Q* = H + 



2h 



'2m 



■d x W, 



i.e. the one with a partner potential 



V' = W 2 + 



'2m 



--d T W 



has Wigner function *-genstates of the same energy as those of H. Specifically, 

H*f = Q**Q*f = f*Q**Q = Ef 
implies that the real functions Q* f *Q* are *-genfunctions of H' with the same eigenvalue E: 

H'*(Q*f*Q*) = Q*Q**Q*f*Q*=E(Q*f*Q*), 
unless f is the Wigner function corresponding to ipo, since Q * fo = 0. 



(43) 
(44) 
(45) 
(46) 



In consequence, E' n = E n+ \ for n > 0. Conversely, for g *-genfunctions of H' , Q* * g ★ Q are *- 
genfunctions of H with the same eigenvalues. 

Moreover, ip' = 1/ipo will be an invalid zero mode eigenfunction of H', as seen from the sign flip in 
( f4l"|) and (44). Consequently, an unnormalized, runaway zero energy solution of the Schrodinger equation 
with V'{x) will invert to the legitimate ground state of H and will permit construction of V given V' . 

For example, starting from the trivial potential with a continuous (unnormalizable) spectrum, 

V = 1, (47) 

and the solution 



% = cosh( — - — ), 



W = tanh(- 



'2mx, 



(48) 



results via ( ^0|) in the symmetric, reflectionless Poschl- Teller potential [15], V = 1 — 2 / cosh 2 ( ^ 
Conversely, starting from this potential, 



V(x) 



1 



cosh 2 (^P)' 



(49) 



there is a single bound state (normalizable to f ip 2 , = 2), 



= sech(- 



W = tanh(- 



'2mx, 



so that 



V' = l. 



Thus, the Wigner function ground state (for m = 1/2) is 



fo(x,p) = — dy 



-iyp 



- dy 



cos(yp) 



2tt J 2 cosh(x/h — y/2) cosh(x/?i + y/2) it Jo cosh(2x/?i) + cosh(y) 

sin(2xp/h) 
sinh(2x/?i) sinh(-7rp) 

(N.B. Not positive definite, nor a function of just H(x,p).) It may be verified directly that 



(50) 
(51) 

(52) 



f tfi x % \ 

Q* fo = (p - y 9x ~ itan H^ + -jdp) j fo(x,p) = 0. 



(53) 
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This appendage of bound states to a potential generalizes [16| to the hierarchy associated with the KdV 
equation. Specifically, 

W(n) = ntanh(^|^) (54) 
connects the reflectionless Poschl- Teller potential 

2 t iw , 2/ \/2mx . , 2 ,2/V / 2mx N 



V'(x)=n — n(n — l)/cosh ( — - — ) to its contiguous V(x) = n —n(n + l)/cosh ( — - — ), (55) 

which has one more bound state (shape- invariance). Recursively then, one may go in TV steps, with the 
suitable shifts of the potential by 2n — 1 in each step, from the constant potential to 

V(N;x) = N 2 - N{N + 1)/ cosh 2 (^^). (56) 

Shifting this potential down by N 2 assigns the energy E = —N 2 to the corresponding ground state 
ipo(N) = sech N (x) (unnormalized) , which is the null-state of ^== d x + W(N). The corresponding (unnor- 
malized) Wigner function is the *-null state of Q(N), 

f (N;x,p) = - [ dy — - . „ N = , 1 ( ,Jt. f (l]x,p), (57) 

7rJo (cosh(2x/ft) + cosh(y)) JV (N — 1)1 \2 smh(2x/n,j / 

where the integral only need be evaluated from the above fo(l;x,p). Alternatively, 

fo(N;x,p) = (sech(x/n)*) 7V - 1 /o(l;x,p) (Wi^/ft))^ 1 . (58) 



The (unnormalized) state above the ground state at E = —(N — l) 2 is yj^dx ~ W(N)j ipo(N — 1), 
and its corresponding Wigner function (setting m = 1/2) is found recursively from the ground state of 
H(N - 1), through Q*(N) * f Q (N - 1) * Q(N), 

p * f Q (N -l)+iN tanh(|) * / (JV - 1)) *Q(N) = (p*f (N-l) + j^j P * fo(N - 1)) *Q(N) (59) 

= (^£y) 2 P*/o(tf-l)*P , 

by virtue of 

g(jv-i)*/o(^-i) = o = / (jv-i)*g*(jv-i). (eo) 

The state above that, at E = —(N — 2) 2 , is found recursively through 

Q*(N)*Q*(N- l)*f (N-2)*Q(N- 1)*Q(N) , (61) 

and so forth. Thus, the entire Wigner *-genfunction spectrum of H(N) is obtained with hardly any reliance 
on Schrddinger eigenfunctions. 



6 Canonical Transformation of the Wigner Function 



For notational simplicity, take fi = 
Canonical Transformations 



1 in this section. The area element in phase space is preserved by 
(x,p)~(X(x,p),P(x,p)) (62) 
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which yield trivial Jacobians {dXdP = dxdp {X, P}) by preserving Poisson Brackets 

_ dudv du dv 
' xp dx dp dp dx 

They thus preserve the "canonical invariants" of their functions: 

{X,P} xp = l hence {x,p} XP = l. (64) 

Equivalently, 

{x,p} = {X,P}, (65) 
in any basis. Motion being a canonical transformation, Hamilton's classical equations of motion are 



preserved, for TC(X,P) = H(x,p), as well [ 17[| . What happens upon quantization? 



Since, in deformation quantization, the Hamiltonian is a c-number function, and so transforms "classi- 
cally", 7i(X,P) = H(x,p), the effects of a canonical transformation on the quantum *-genvalue equation 
of Lemma 1 will be carried by a suitably transformed Wigner function. Predictably, the answer can be 
deduced from Dirac's quantum transformation theory. Consider the canonical transformations generated 
by F(x,X): 



Following Dirac's celebrated exponentiation |18[ of such a generator, in the implementation of |12| , [191] , 
the energy eigenfunctions transform canonically through a generalization of the "representation-changing" 
Fourier transform. Namely, 

M*) =N E jdX e lF ^ V E (X) . (67) 



Thus, 

i 2 



f(x } p) = Jdy J dX l e -iF*(*-v/W e~ iyp J dX 2 e iF (*+y/W y E (x 2 ) . (68) 

The pair of Wigner functions in the respective canonical variables, f(x,p) and 

F(X, P) = ±-JdYt>*{X- ^Y) e-* YP *(X + \y), (69) 
are connected by a transformation functional %(x,p; X, P), 

f(x,p) = JdX JdPZ(x,p;X,P) ®F(X,P)= JdX fdP%(x,p;X,P) F(X,P) , (70) 

where © is with respect to the variables X and P. 

To find this functional, let X = \{X X + X 2 ) and Y = X 2 - X x , so that J dXi J dX 2 = J dX JdY. 
Noting that 

**(X - |y) V(X + \y) = J dP e lYP F(X, P) , (71) 
it follows that (|6^) reduces to 

f( Xj p) = L^i! J dy fdX 1 e- iF *( x - y / 2 > Xl H*(X 1 )e- iyp JdX 2 e iFix+y/2 ' X2) y{X 2 ) (72) 



\Nl 
2vr 



JdXdYdy e~ iyp e -« ?, (*-i//2.^-y/2) _ y/2) ^{X + Y/2) ^+v/2,x+Y/2) 

dXdPdYdy e -iyP+iPY-iF*(x-y/2,X-Y/2) + iF( X +y/2,X+Y/2) jr( X ,P) , 
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which leads to 

Lemma 3: %(x,p; X, P) = WlfdYdy exp (-iyp + iPY - iF*(x - §, X - |) + iF(x + f , X + |) 

Corollary 2: T/iis phase-space transformation functional obeys the "two-star" equation, 

H(x,p)*Z(x,p;X,P) =%{x,p;X,P) ®H{X,P) , (73) 

as follows from i?(x, —id x ) exp(zF(x, X)) = W(X, idx) exp(iF(x, X)). If J 7 satisfies a @— genvalue equa- 
tion, then / satisfies a *-genvalue equation with the same eigenvalue, and vice versa. ■ 

Note that, by virtue of the spectral projection feature (1€ Jl9|), this equation is also solved by any 
representation-changing equal-energy bilinear in real Wigner *-genfunctions of H and 7i, 

<I(x,p;X,P)=Y / 9(E) f E (x,p) F E (X,P) , (74) 

E 

for arbitrary real g(E). Such a bilinear transformation functional is nonsingular (invertible) if and only if 
g(E) has no zeros on the spectrum of either Hamiltoniani. 

As an example, consider the linear potential again, which transforms to a free particle (TC = P) through 

F = --X 3 -xX =^ P = ~X, x = P-X 2 . (75) 
3 

By direct computation, 

%(x,p;X,P) = 2 2 ^Ai{2 2 '\x + X 2 -P)) 6{p + X) = (2n) 2 J dE f E (x,p)T E {X,P) S(p + X) . (76) 

Note Ne = \j\f2m for the free particle energy eigenfunction normalization choice ^ E (X) = (2tt)~ 1 ^ 2 exp(iEX). 
Thus, indeed, the free particle Wigner function, Te{X,P) = 5(E — P)/(2ir), transforms to 

1 r 9 2 / 3 / \ 

f(x,p) = —J dPdX % 5{E -P) = —Ai [2 2 l\x +p 2 - E)) , (77) 



as it should; and (73) is seen to be satisfied directly by virtue of the linearity of the respective Hamiltonians 



in the variables P, x, conjugate to those of the arguments of 5(p + X). 

The structure of the result in ([76]) underscores that the linear potential is as "close to classical" as one 



can get, in simple quantum mechanics. It has been noted before [12| that the transformation functional 
for linear potential wave functions is exactly the exponential of the classical generating function for the 
canonical transformation to a free particle, and that this is not the case for any other potential. The 
present result for the transformation functional for Wigner functions is further evidence for this "close to 
classical" behavior. The delta function 5(p + X) in ([?(]) is half 'of the classical story. Were the Airy function 



2 In general, if the transformation functional effects a map to a free particle, the P integration is trivial in (^), and the 
result for the Wigner function of the x,p theory is just an average over X of the transformation functional. That is, if 
J-{X, P) — S(P — k(E)), where k(E) is the momentum-energy relation for the free particle theory in question, 

f{x,p) = J dX JdP %(x,p;X,P) F(X,P) = JdX %(x,p;X,k(E)). 

One might then be tempted to wonder if just %(x,p;X, P) = ip* P (x — hX/2) e~ zXp tpp(x + hX/2)/2-K = (3(x,p; X, P). However, 
what determines the allowed range for PI It is always possible to embed any real energy spectrum into the real line, but 
knowing this does not help at all to determine what points are to be embedded. From the point of view of this paper, even 
when the spectrum is obvious, such a choice for the transformation functional in general does not satisfy the two-* equation 
(|7^). Rather, the equation fails by total derivatives that vary contingent on particularities of the case. E.g., for free-particle 
plane waves, iPe{x) = exp(iEx), so that p*(J5 — (J5@P = dx®- This choice for %, then, does not yield useful information 
on the Wigner functions. 
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also a delta function of its argument, we would have an exact implementation of the X, P \— > x, p classical 
correspondence. As it is, there is some typically quantum mechanical spread around the classical constraint 
x + X 2 — P = 0, in the form of oscillations of the Airy function, and, in consequence, the Wigner functions 
of the free particle do not retain their delta-function form under the canonical transformation to the linear 
potential Wigner functions. Reinstating h into (p6|)i, and taking the limit fi — > converts the Airy function 
to a delta function, 5(x + X 2 — P), thereupon producing the complete classical correspondence between 
the two sets of phase space variables, in that limit. 

As already seen, there is substantial non-uniqueness in the choice of transformation functional. For 
example, for the linear potential again, (|73|), 

(x+p 2 )*6(x,p;X,P)) = &(x,p;X,P)) ®P (78) 

is also satisfied by a different (and somewhat simpler) choice: 

&(x,p;X,P) =exp(-i(^X 3 + 2(x+p 2 -P)X)) . (79) 

This transformation functional also converts the free particle Wigner function, Te{X, P) = 5{E — P)/2tt, 
into an Airy function (as above) after integrating over the free particle phase space, JdXdP. 

Actually, it is not necessary to integrate over the phase space. In general, ^-multiplying a delta function 
spreads it out, and yields a Fourier transform with respect to the conjugate variable. Thus, for the example 
considered, 

eJ (-|X 3 -2(xV-PW i(5 (p_ £ ) _ e 2iX(P-E)}_ I dZe -2iZ(P-E) e i(-lz3-2{x+p*-P)Z) 

vr J 

= e 2iX(P-E)l f dZe i(-lz3-2( X+p i-E)Z) = e 2iX(P-E) 2 2/3 Ai Q2/3( x+p 2_ E)) (8Q) 
vr J 

Hence, 

JdX [dPe i (-i x3 - 2 ( x+ P 2 -^ x U5(P-E)= 2 2 /\ Ai(2 2 / 3 (x+p 2 -E)) . (81) 
Compare this to the action of the above %(x,p; X, P), 

(Ai(2 2 / 3 (x + X 2 -P))5(p + A)) *6(P-E) = (82) 
e 2iX ( p ~ E ^- [dZe- 2iZ ( p - E Ui(2 2 / 3 (x + Z 2 -P))5(p + Z) = e n P +x)(P-E) I Ai( tf/3, x + p 2 _ p)) 

IT J IT 

Aside from innocuous normalizations, the difference in the two transformation functionals acting on the 
free particle Wigner function is just the phase factor e 2vp ^ p ~ E \ and the argument of the Airy function, 
where E has been replaced by P. Indeed, the phase factor precisely compensates for the different energy- 
eigenvalue occurring in the argument of Ai, when acted upon by (x + p 2 )*. Such simple phase factors may 
be used to shift a *-genvalue whenever the Hamiltonian is linear in any variable. 



7 Illustrations using Liouville Quantum Mechanics 

Summary illustration of all the above, in particular the canonical transformation effects on Wigner func- 
tions, is provided by the Liouville model |20[j. Our conventions for the model (which are essentially those 



of [21], with their m = l/(4vr) and their g = 1) are given by 

Huouville = p 2 + e 2x . (83) 



3 The exponent of the integrand turns to iy(E — x — p 2 — h 2 y 2 /12). 
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The energy eigenfunctions are then solutions of 

/ d 2 



v -^ + e2X j^(x)=Eip E (x) . (84) 
The solutions are Kelvin (modified Bessel) K functions, for < E < oo, 

Mx) = ~ ^/sinh (ttVe) K is/ z(e x ) , (85) 

which are normalized such that dx i/j e (x) ipE 2 ( x ) = S(Ei — E%). There is no solution [2C] for E = 0. 

For completeness, consider the Fourier transform (including a convergence factor, necessary for x — > — oo 
to control plane wave behavior, but not for x — > oo): 

/+oo 
dxe- tx{p+le) ^ E (x) (86) 
-oo 

- JsmUnVE) 2~*»*> T ( zlk+lll+l^ ) T ( -i (p + ie) - _ 



^ /sinh^E) 2- ^ 2 ry 2 

This follows, e.g., from a result in p^], Vol II, p 51, eqn (27): 

+ ~ dz ^„ W = 2,- r (1±|±^) r (i±|^) , (87, 

valid for 3?(1 + ± v) > (i.e. the previous transform is valid for e > 0). The right hand side of this 
last relation clearly displays the symmetry v — ► — u, which just amounts to the physical statement that 
the energy eigenfunctions are non-degenerate for the transmissionless exponential potential of the Liouville 
model. 

Further note the effect on §e{p + ie) of shifting p — > p + 2i, using T(l + z) = z T(z), 

. f-i(p + ie) +iVE\ f-i(p + ie) - %Ve\ , , 
<f> E (p + 2i + ie) = 4(^^ J- ^— j ( — HI J- ^U E (p + ie ) 

= (E-{p + ief) $ E {p + ie) ■ (88) 

So, as e — > 0, $e(p + 2i) = (E — p 2 ) $e(p). But this simple difference equation is just the Liouville energy 
eigenvalue equation in the momentum basis, 

(p 2 - E) $ £ (p) + e 2id " § E {p) = . (89) 

Such first order difference equations invariably lead to gamma functions [23]. Below, it turns out that the 
Wigner functions also satisfy momentum difference equations, but of second order. 

Many, if not all, properties of the Liouville wave functions may be understood from the following integral 
representation |24|| , Ch VI, §6.22, eqn (10). Explicitly emphasizing the abovementioned non-degeneracy, 

1 r+oo 

K ik {e x ) = K_ ik {e x ) = - e nk/2 / dXe ieXsinhX e ikX . (90) 

2 J — OO 



(Also see [25], 9.6.22.) This integral representation may be effectively regarded as the canonical trans- 
formation of a free particle energy eigenfunction, e tkX , through use of the generating function F(x,X) = 
e^sinhV. Classically, p = dF/dx = e^sinhV, and P = —dF/dX = — e^coshX, so P 2 — p 2 = e 2x . That 
is, HLiouviiie = ~Hfree = P 2 under the classical effects of the canonical transformation. The quantum effects 
are detailed below, by *-acting with the Liouville and free Hamiltonians on the suitable transformation 
functional. 
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The Liouville Wigner function may be obtained from the definition ([!]) in terms of known higher 
transcendental functions, 



f(x,p) 



1 



I f+oo 

2^ 



dy ^ smh(WE) K t ^{e x -y' 2 ) e'™ K^e^' 2 ) (91) 
1 + 2i y/E 1 - 2iy/E 1 + 2iy/E + 4ip 1 - 2iy/E + Aip \ 



The following K-transform was utilized to express this result in closed form, 



dw (wz) 1/2 vf" 1 K„ {a/w) K v (wz) = 2- ,J - 5 /V G\ 



10 
01 



a 2 z 2 



The right hand side involves a special case of Meijer's G-function 



p — a 



-fj, — cr 1 v 1 f 
"2 '4 + 2 '4 ~ 2 



dj, i=l,...,p 
bj, j = l,...,q 



(92) 
(93) 

(cf. p2fl , §5.3), which is fully symmetric in the parameter subsets {ai, . . . , a n }, {a n+ i, . . . , a p }, {hi, . . . , 6 m }, 
and {6 m+ i, . . . , It is possible to re-express the result as a linear combination of generalized hyperge- 
ometric functions of type 0F3, but there is little reason to do so here. This transform is valid for 3fta > 0, 
and is taken from [26], p 711, eqn (55).Q The transform is complementary to [27], §10.3, eqn (49), in an 
obvious way, a K-transform which appears in perturbative computations of certain Liouville correlation 
functions lElfl. 



The result ( |9~l"D may be written in slightly different alternate forms, 

1 + 2i\[E -2ip 1 - 2i\/E - 2ip 1 + 2i\[E + 2ip 1 - 2ixfE + 2ip 



smh(ir\/E)e x 40 / e 4x 
-G, 



47T 3 



T 4 



16 



sinh(7rv E) 
8^3 



r 4 



16 



i\[E — ip —i\[E — ip i\[E + ip —i\jE~-Y 



ip 



(94) 



by making use of the parameter translation identity for the G-function (p2[ §5.3.1, eqn (8)): 



z x G™ z 



a r 



U pq Z 



a r + A 
b s + X 



(95) 



Yet another way to express the result utilizes the Fourier transform of the wave function, (|8q), in terms 
of which the Wigner function reads, in general, 



f(x,p) -- 

The specific result ( |S6"1 ) then gives, as e — » 
f( x ,p) = ( — \ sinh(vrV^) / dk e ixk 4~ i( - k/2 



1 \ 2 r +co i t. 1 

—J y ^ dfe^(p--A;)e- fc ^( p+ _fc). 



+fe) r ( j (P ~ k / 2 ~ ie ) ~ 



(96) 



(97) 



/ i (p - k/2 - ie) + / -j (p + fc/2 + ie) + / -i {p + fc/2 + ie) - v/e \ 



However, this is a contour integral representation of the particular G-function given above. Because of 
the e prescription, the contour in the variable z = k/2 + ie runs parallel to the real axis, but slightly 



4 There is an error in this result as it appears in EtJ], Vol II, §10.3, eqn (58), where the formula has a 2 z 2 /4 instead of 
a 2 z 2 /16 as the argument of the G-function. The latter argument is correct, and appears in Meijer's original paper cited here. 
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above the poles of the T-functions located on the real axis at z = p — y/E, z = p + VE, z = —p + \/E, 
and z = —p — y/~E . Changing variables to s = \iz yields 

f(x,p) = sinh(vr^) -L / ds(^r) T ( ip - ^ - s) x (98) 



8vr 3 v 2-Ki Jc \ 16 J \ 2 

XT [-^- S ) T [ 2 S ) T { 2 S ) ' 

where the contour C in the s-plane runs from —zoo to +ioo, just to the left of the four poles on the 
imaginary s axis at i(p + y/E)/2, i(p — yE)/2, i(—p + vE)/2, and i(—p — VE)/2. This is recognized 
as the Mellin-Barnes type integral definition of the G^-function (cf. |22|], §5.3, eqn (1)) in agreement with 
the second result above, (94). 

The translation identity ( |95| ) is seen to hold by virtue of (|98|), through simply shifting the variable 
of integration, s. Moreover, deforming the contour in (|98|) to enclose the four sequences of poles s n = 
n + i(±p =t s/E) /2 reveals the equivalence of this particular G-function to a linear combination of four qF^ 
functions, one for each of the sequences of poles. Evaluating the integral by the method of residues for all 
these poles produces the standard hypergeometric series. 

It should now be straightforward to directly check that the explicit result for f(x,p) is indeed a solution 
to the Liouville *-genvalue equation, 

H Li(mvaie *f(x,p) = ((p - l -d x f + e 2 (-+5^ f^ p ) = E f(x,p) . (99) 

For real E and real f(x,p), the imaginary part of this *-genvalue equation is 

(-pd x + e 2x sin d p ) f(x,p) = , (100) 

while the real part is 

p 2 - E - ~dl + e 2x cos df}j f(x,p) = . (101) 
The first of these is a first-order differential/difference equation relating the x and p dependence, 

e- 2x d x f(x,p) = -L (f(x,p + i) - f(x,p- i)) . (102) 
zip 

Similarly, the real part of the *-genvalue equation is a second-order differential/difference equation, 

e- 2x (p 2 -E--^d 2 )j f( x ,p) + ^f( x ,p + i) + f( x ,p-i)) = 0. (103) 

The previous first-order equation may now be substituted (twice) into this last second-order equation, 
to convert it from a differential/difference equation into a second-order difference only equation in the 
momentum variable, with non-constant coefficients. That is, 

P 2 -E) f( x ,p)+(il\ (f(x,p + 2i)-2f(x,p) + f(x,p-2i)) 

e 2x e 2x 
+i (f(x,P + i) - f(x,p- i)) + — (f(x,p + i) + f(x,p - i)) . (104) 

We leave it as an exercise for the reader to exploit the recursive properties of the Meijer G-function and 
show that this difference equation is indeed obeyed by the result (|91~|). Rather than pursue this in detail, 
we turn our attention to the transformation functional which connects the above result for / to a free 
particle Wigner function. 
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Given (|9C|), it follows that 

i(>e(x) = - \/sinh(7r\/l) A. /p(e x ) = ±- x/sinh(7r\/^) e^ /2 / + °° e ^ sinh * e *^x ( 105 ) 

hence = (4vr v A Ee 7rv/ ^ sinh(vry^)) 1 /2/27r, if we choose a 5{E\ — E2) normalization for the free particle 
plane waves as well as for the Liouville eigenfunctions. Therefore, Lemma 3 yields 

1(x, p; X, P) = I^L / dYdy exp (-iyp + iPY - iF* (x - y/2, X - Y/2) + iF(x + y/2, X + Y/2)) (106) 
2ir J 

= — ^ (inVEe^ sinh(7rv A E)) J dYdy exp {-iyp + iPY - ie x ' y/2 sinh (x - — ^ + ze x ' +?//2 sinh + 

= J_ (^v^e^ S inh(7rv^) ) J d (^^) exp (i (P - p) + sinh (*±I 

X ( d ( Y ~ V \ exp f i (P + p) ^ ~ V + ie x_x sinh ^ ' " 



Y 



2 

We thus conclude, 



%(x,p;X,P) = ^VEe^ S mh(nVE) e^ p K l(P ^ p) (e x+x ) K l(P+p) (e x ~ x ) . (107) 



We now check that this result obeys (73) and, in so doing, carry out the nontrivial steps needed to 
show the Liouville Wigner functions satisfy the Liouville *-genvalue equation (|99|). That is to say, we shall 
show 

Up - \~dx? + e^+^J %(x,p; X, P) = S(x,p; X, P) ((P + \^ xf) , (108) 
or, equivalently, 

({p ~ I'd,? + e 2 ^+t^) - (P + \~dx)^ K l(P _ p) {e x+x ) K l(P+p) {e x ~ x ) = . (109) 

Specifically, 

=1 - ~9x) K i{P _ p) (e* +x ) K l{P+p) (e x - x ) = -e 2x K' i{p _ p) (e* +X ) K' t{p+p) (e^ x ) , (110) 



-ip d x -iPd x j K i{P _ p) (e x+x ) K l(P+p) (e x ~ x ) (111) 
= ~i (p + P) e X+X K' i{ P- P )^ X+X ) K l{P+p) {e*~ x ) -i(p-P) e*~ x K l{P _ p) (e x+x ) K' i{P+p) {e^ x ) , 

and 

e 2 ^ x+ i^^ K i(P _ p) ^+ x )K i{P+p) (e x - x ) = e 2x K 1+l{P _ p) {e j:+x )K. l+l{P+p) {e x - x ) . (112) 

Now, recall the recurrence relations ( p5| , 9.6.26) 

K 1+l[P _ p) {e* +x ) = -K' l{P _ p) (e x+x ) + i (P-p)e- x - x K l{P _ p) (e x+x ), (113) 

K_ 1+l(P+p) (e*- x ) = -K' i(P+p) (e x - x )-i(P + p)e- x+x K iiP+p) (e x - x ). (114) 
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So the previous relation (112) becomes 

e 2(,+f ~d P ) Kl(p „ p) (e* +x )K l{P+p) (e*- x ) = e 2x K' i(P _ p) (e x+x ) K' l{P+p) {e x ~ x ) (115) 

+i (P + p) e* +x K' i{P _ p) (e* +x ) K i(P+p) (e x - x ) 
-i (P - p) e*~ x K l(P _ p) (e x+x ) K' {P+p) (e x - x ) 
+ (p 2 - P 2 ) K l(P _ p) (e x+x )K l{P+p) (e x - x ). 

The sum of ( |110| ), ( |111| ), and ( 115 ), shows that ( |109| ) is, indeed, satisfied. 

Integrating over X and P the product of %(x,p;X,P) and the free particle Wigner function, as given 
here by (4ir\/ r E)~ 1 5(P — \^E), yields another expression for the Liouville Wigner function which checks 
against the previous result, (^). Using (92) and the parameter translation identity for the G- function, 
this other expression is just (| 



Supersymmetric Liouville quantum mechanics is obtained by carrying through the Darboux construc- 
tion detailed above (with h = 1 = 2m), for the choice 

W{x) = e x . (116) 

Conventions essentially follow |2S[| . 

The first Hamiltonian of the essentially isospectral pair is then 

H = p 2 + e 2x — e x , (117) 

and the allowed spectrum is < E < oo , including zero-energy, for which there is a bounded wave function 
normalized as part of the continuum, 

The other, E > 0, eigenfunctions are 

M*) = (^f e * cosh 0^)) V2 fa-ivst*) + K frvs(*)) > (H9) 

again normalized so that dxt/; El (x) tpE 2 (^) = S(E\ — E2). 
The second Hamiltonian of the pair is 

H' = p 2 + e 2x + e x , (120) 
and the allowed spectrum is < E < 00, excluding zero energy!. The E > eigenfunctions are then 

M*) = (l^ eXcosh (^)) 1/2 ( iK h-iVB^)-iK h+iVB (e^)) , (121) 

and may be obtained from the previous E > eigenfunctions, as ip' E (x) = (d x + W) iJ)e{x)- 

For both Hamiltonians, the Wigner functions are straightforward to construct directly, once again 
leading to the K-transform (^) and particular Meijer G-functions. We find it sufficient here to consider 
only the ground state for H, 

1 z'+oo 9 

M*,P) = ^2 J ^ dye- 2eXcosh ^-^ = -K 2ip (2e x ) , (122) 



3 The candidate ip' (x) = l/ipo(x) — y^n exp(e x ) solves the Schrodinger equation, but is obviously unbounded, as expected. 
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a single modified Bessel function. It smoothly satisfies (p — iW(x)) * /o = 0, and hence the *-genvalue 
equation H * f = 0. 
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